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ABSTRACT 


i 
l 
I 

V 

An  experiment  to  determine  the  interaction  of  an  intense 
electromagnetic  pulse  (EMP)  ,  such  as  that  produced  by  a 
nuclear  detonation  above  the  earth's  atmosphere,  was  performed 
in  March,  1986  at  Kirtland  Air  Force  Base  near  Albuquerque, 

New  Mexico.  The  results  of  that  experiment  have  been 
publishecy^J^anT without  analysis.  Following  an  introduction 
of  the  corona  phenomenon,  the  reason  for  interest  in  it,  and  a 
review  of  the  experiment,  this  paper  discusses  five  different 
analytic  corona  models  that  may  model  corona  formation  on  a 
conducting  line  subjected  to  EMP.  The  results  predicted  by 
these  models  are  compared  with  measured  data  acquired  during 
the  experiment  to  determine  the  strengths  and  weaknesses  of 
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EXECUTIVE  SUMMARY 


A  nuclear  detonation  in  or  above  the  Earth's  atmosphere 
produces  an  intense  electromagnetic  pulse  (EMP) .  The  inter¬ 
action  of  this  EMP  with  the  nation's  power  systems  could 
disrupt  these  systems  leaving  large  areas  of  the  country 
without  electric  pcwer.  This  is  a  problem  of  present  concern  [1], 
The  extent  and  duration  of  these  system  blackouts,  should  they 
occur,  is  presently  unknown.  To  determine  the  probability 
of  such  an  occurrence  and  its  severity,  the  Department  of 
Energy  (DOE)  has  formulated  a  program  through  the  Oak  Ridge 
National  Laboratory  (ORNL)  to  assess  the  effects  of  EMP  on 
commercial  power  systems. 

The  major  research  effort  in  this  program  is  being 
conducted  by  the  Westinghouse  Electric  Corporation.  This 
effort  has  produced  several  computer  models  for  EMP  interaction, 
has  performed  experiments  on  power  systems  eq  ipment  to  charac¬ 
terize  its  EflF  behavior,  and  has  developed  an  assessment 
methodology.  The  methodology  uses  fault  tree  analysis  to 
estimate  system  behavior. 

In  the  area  of  numerical  modeling  of  the  EMP  interaction 
of  power  lines  with  the  EMP,  a  number  of  models  have  been 
developed.  These  models  are  generally  linear  models  and  the 
effects  of  corona;  i.e.,  ionization  of  the  air  in  the  vicinity 
of  the  conductor,  are  not  taken  into  account.  Corona  is  known 
to  alter  the  shape  and  magnitude  of  an  injected  pulse  by  a 
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damping  action.  However,  the  models  used  to  predict  corona 
effects  on  injected  pulses  are  not  appropriate  for  modeling 
an  EMP  coupled  pulse.  These  injected  pulse  corona  models 
generally  use  a  pulse  ri3e  time  of  microseconds  as  compared 
to  a  few  tens  of  nanoseconds  rise  time  for  the  EMP  coupled 
pulse.  Also,  the  EMP  pulse  is  characterized  as  a  distributed 
source  which  requires  different  modeling  techniques.  There 
was  no  known  data  on  corona  effects  for  fast  rising  (nanosecond) 
impulse  and  distributed  excitation.  Consequently,  in  March, 

1986,  an  EM?  corona  experiment  was  conducted  at  Kirtland  Air 
Force  Base,  New  Mexico.  This  work  was  funded  in  part  by  the 
DOE  and  by  the  Defense  Nuclear  Agency.  This  report  examines 
some  of  the  experimental  results  in  detail;  viz,  a  baseline 
test  case  representative  of  the  data  collected,  and  compares 
these  results  with  five  corona  models. 

In  general,  three  of  the  models  fit  the  measured  data 
reasonably  well.  In  all  cases,  the  limitations  of  the  experi¬ 
mental  configuration  presented  some  difficulty  in  analyzing 
the  data.  However,  these  difficulties  were  expected  and 
could  not  be  avoided  using  the  present  facilities.  In  particular, 
the  experimental  configuration  allows  the  examination  of  only 
"positive"  conductor  charge  and  has  a  limited  "clear  time"  of 
about  90  ns.  "Clear  time"  is  the  time  prior  to  the  appearance 
of  significant  reflections  in  the  experimental  configuration. 
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The  results  of  the  tests  indicate  that  corona  is  present 
under  BMP  conditions.  In  fact,  the  experimental  team  was 
unable  to  find  a  configuration  without  corona.  Therefore, 
all  comparisons  to  the  no  corona  case  are  with  respect  to  a 
calculated  no  corona  response.  While  all  models  fit  the 
data  over  certain  ranges,  no  single  model  is  capable  of  covering 
the  entire  data  range.  It  is  clear  that  ignoring  corona  in 
the  evaluation  of  EMP  impact  will  produce  a  conservative  result; 
however,  the  degree  of  conservatism  remains  in  question. 
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I .  INTRODUCTION 


The  interaction  of  an  intense  electromagnetic  pulse 
(KMP) ,  such  as  that  produced  by  a  nuclear  detonation  in  or 
above  the  earth's  atmosphere,  with  a  power  system  is  a  problem 
of  present  concern  [1].  The  EXP  could  disrupt  the  functioning 
of  the  system  and  leave  large  areas  of  the  country  without 
electrical  power.  Although  occasional  localized  power 
disruptions  are  experienced  in  commercial  power  systems,  the 
large  geographical  coverage  of  a  high  altitude  EMP  suggests 
that  there  is  a  possibility  that  power  outages  over  extended 
areas  might  occur  in  the  event  of  such  a  nuclear  detonation. 

The  extent  of  such  a  power  black-out  (should  it  occur) , 
{.s  well  as  the  length  of  time  that  the  black-out  might  exist, 
is  presently  unknown,  in  order  to  determine  the  probability 
of  such  an  occurrence  and  its  severity,  should  it  occur,  the 
U.  S.  Department  of  Energy  (DOE)  has  instigated  a  program 
through  the  Oak  Ridge  National  Laboratory  (ORNL)  to  assess  the 
effects  of  an  EMP  on  the  commercial  power  system.  This 
effort,  being  conducted  by  the  Westinghouse  Electric  Co.,  has 
developed  a  number  of  computer  models  for  EXP  interaction 
studies,  has  performed  experiments  on  power  system  equipment 
to  characterize  their  EMP  behavior,  and  has  developed  an 
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assessment  methodology  based  on  a  fault  tree  analysis  to 
estimate  the  system  behavior. 

In  the  area  of  numerical  modeling  of  the  EMP  interaction 
of  power  lines  with  the  EMP,  a  number  of  transmission  line 
coupling  models  have  been  developed.  The  majority  of  these 
are  linear  models,  and  consequently,  the  effects  of  corona 
(either  due  to  the  energized  line,  or  to  the  EMP  itself)  are 
not  taken  into  account.  However,  it  is  desirable  to 
understand  the  effect  that  corona  has  on  the  coupled  EMP 
response  of  long  lines  and  to  ascertain  if  completely 
neglecting  the  corona  effects  in  the  power  system  assessment 
is  reasonable. 

Several  investigators  have  previously  studied  the 
effects  of  corona  on  the  EMP  response  of  long  lines 
analytically.  Much  of  this  work,  however,  has  not  been 
directly  verified  by  experimental  observation.  Consequently, 
in  March,  1986,  an  EMP  corona  experiment,  funded  in  part  by 
DOE  through  ORNL  and  by  the  Defense  Nuclear  Agency  (DNA) ,  was 
conducted  at  Kirtland  Air  Force  Base  in  New  Mexico.  The  data 
obtained  in  this  measurement  program  have  been  published 
[2,  3],  without  analysis  comparing  the  experimental  results 
with  the  various  corona  models.  The  purpose  of  the  present 
document,  therefore,  is  to  examine  some  of  the  experimental 
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results  in  more  detail  and  to  compare  them  with  the  results 
from  the  calculational  corona  models. 

In  Section  II  of  this  report,  brief  discussions  of  the 
formation  and  the  effects  of  corona  are  presented.  This  is 
followed  in  Section  III  by  a  review  of  the  experiment  and  the 
baseline  test  configuration  chosen  for  comparing  the 
predictions  of  several  corona  models  with  the  measured  data. 
The  rationale  for  selecting  the  baseline  test  configuration  is 
also  summarized. 

After  these  general  background  sections,  a  description 
of  the  different  corona  models  analyzed  in  this  paper  is  given 
in  Section  IV.  The  results  predicted  by  each  of  the  corona 
models  are  compared  with  the  data  acquired  for  a  chosen 
baseline  configuration  and  an  determination  of  the  appropriate 
parameters  for  each  model  is  made  to  determine  the  best  fit  of 
the  model  to  the  test  case  data. 

Following  this,  Section  V  concludes  the  report  with  a 
summary  of  the  important  observations  and  some  possible 
directions  for  future  research. 
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II.  WHAT  IS  CORONA?  *  A  BRIEF  BACKGROUND 

In  this  section,  an  overview  of  Peek's  law,  describing 
the  electric  field  limit  beyond  which  corona  formation  will 
begin,  is  presented.  This  is  followed  by  a  brief  discussion 
of  the  effects  of  corona  on  power  lines. 

A .  Peek 1 s  Law 

Typically,  the  atmosphere  is  a  dielectric,  or  insulator, 
for  which  the  nominal  dielectric  strength  is  approximately  30 
kV/cm.  If  the  local  electric  field  exceeds  this  value,  corona 
will  begin  to  form.  The  actual  value  of  the  dielectric 
strength  may  change,  however,  depending  upon  wire  roughness, 
barometric  pressure  and  temperature,  as  indicated  by  Peek's 
law  [4,  5]: 


Ec  =  E0  m  6  t1  +  V(5r)1/2] 


for  coaxial  geometry 


where : 


D  is  the  threshold  electric  field  above  which 

C 

corona  may  be  expected  to  form  . 
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Eq  is  the  sparking  field  strength  for  air  under 

STP  and  uniform  field  for  gaps  on  the  order 
of  1  cm.  An  approximate  value  is  31  kV/cm. 

m  is  a  coefficient  describing  the  conductor 
surface  roughness:  0.6  <  m  <  1.0  ,  with 

m  =  1.0  for  smooth,  round  conductors 

r  is  the  conductor  radius  (inner  conductor  for 
a  coaxial  geometry) . 

For  r  in  cm,  the  value  of  k  is  0.308  cm1^2  .  The 
variable  6  is  the  air  density  factor  relative  to  that  at 
standard  temperature  and  pressure  (STP) ,  which  is  equal  to 
0 . 392p/T  ,  where  p  is  the  pressure  corrected  to  sea  level  in 
mm-Hg  and  T  is  the  ambient  temperature  in  degrees  Kelvin. 
This  law  is  valid  over  a  large  range  of  pressures,  for 
0.1mm  <  r  <  several  cm,  and  for  applied  voltages  from  DC 

to  several  kHz. 

The  ambient  air  pressure,  corrected  to  sea  level,  and 
temperature  data  on  the  date  of  the  baseline  test  were 
obtained  from  the  National  Weather  Service  office  in 
Albuquerque,  NM.  The  nominal  values  for  these  parameters  for 
the  baseline  case  chosen  are  as  follows: 

Eq  =  3.1  x  106  (Volts/meter) 

0.6  <  m  <  1.0 

752.1  <  p  <  754.9  (mm-Hg) 
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279.1  <  T  <  280.2  (degrees  Kelvin) 
r  =  1.07  (cm) . 

Figure  1  contains  a  plot  of  Ec  as  a  function  of  m 
for  the  ranges  of  p  and  T  indicated  above. 

B.  Corona  Effects 

When  the  local  electric  field  intensity  exceeds  the 
dielectric  strength  of  the  atmosphere,  ionization  of  the  air 
occurs  in  the  region  of  high  field  intensity.  This 
effectively  transforms  the  air  in  that  region  into  a  weak 
electrical  conductor.  It  is  well  known  that  one  way  the  local 
electric  field  may  be  enhanced  is  by  the  presence  of  a  wire 
with  a  relatively  small  radius  carrying  a  large  current.  When 
the  air  is  ionized  around  a  wire  because  of  the  enhanced  local 
electric  field,  corona  (from  the  Greek  kronos,  or  crown)  is 
said  to  have  formed.  The  formation  of  corona  on  a 
transmission  line  is  generally  undesirable,  since  it 
dissipates  some  of  the  line's  power  and  changes  the 
capacitance  of  the  transmission  line,  suppressing  the  amount 
of  current  the  line  may  carry. 

Studies  of  the  coupling  of  an  electromagnetic  pulse  (EMP) 
produced  by  a  nuclear  detonation  to  a  lohg  wire  indicate  that 
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Figure  1.  Breakdown  Electric  Field  Strength  as 
Predicted  by  Peek’s  Law  [4,5]  for 
0.5  <  m  <  1.0  on  Baseline  Test  Date 


large  electrical  currents  may  be  induced  on  the  wire.  These 
currents,  in  fact,  may  be  large  enough  to  generate  corona  on 
the  line,  in  much  the  same  way  large  currents  on  power 
transmission  and  distribution  lines  do.  However,  because  of 
the  transient  nature  of  EMP,  the  EMP  induced  corona  is  also 
short  lived.  Nevertheless,  it  may  still  reduce  the  amount  of 
current  induced  on  the  line  by  the  EMP. 

Before  presenting  the  details  of  this  analysis,  a 
description  of  an  experiment  performed  to  measure  EMP  induced 
corona  effects  is  given  in  Section  III.  Several  analytical 
models  of  the  effects  of  EMP  induced  corona  on  the  induced 
line  current  have  been  developed.  These  are  discussed  in 
Section  IV.  Following  each  discussion,  the  predictions  of  the 
analytical  model  are  compared  with  experimental  results.  A 
separate  sub-section  then  compares  each  model  in  terms  of  how 
well  its  predictions  fit  the  experimental  data.  Finally, 
Section  V  presents  concluding  remarks  and  describes  possible 
directions  for  future  research. 
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III.  DESCRIPTION  OF  EXPERIMENT  AND  BASELINE  TEST  CASE 


In  this  section,  a  general  overview  of  the  experiment  and 
the  processing  of  the  data,  following  its  acquisition,  are 
outlined.  The  baseline  test  case  is  described  in  this 
context . 

A.  Experimental  Configuration 

A  pre-test  survey  of  the  EMP  simulators  available  to 
conduct  this  experiment  indicated  the  VPD-II  at  Kirtland  Air 
Force  Base  in  Albuquerque,  New  Mexico  as  the  simulator  with 
the  most  appropriate  field  and  working  volume.  Since  a  common 
variable  in  at  least  three  corona  models  (Baum  [6],  Townsend 
[7,  8],  and  conductivity  [7]  models)  studied  prior  to 
performing  this  experiment  was  the  angle  of  E-field  incidence 
with  respect  to  the  conducting  line,  this  experiment  consisted 
of  a  test  line  suspended  at  various  known  angles  with  respect 
to  the  vertical  electric  field  emitted  by  the  VPD-II.  Induced 
currents  were  measured  with  an  in-line  current  derivative 
(dl/dt)  sensor,  designed  by  EG&G  and  designated  an  OMM-2  probe 
[93- 


Other  variables  explored  during  this  experiment 
included:  test  wire  length,  effective  wire  diameter,  changes 
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in  surface  roughness,  the  addition  of  a  dielectric  sheath 
(insulation)  to  the  test  line,  and  variations  in  field 
strength  of  the  pulse  emitted  from  the  VPD-II. 

Figure  2  indicates  the  general  configuration  used  in  this 
test.  A  conducting  line  was  suspended  from  VPD-II  support 
pole  #8  at  a  chosen  height,  H,  and  extended  diagonally  over 
the  ground  plane  to  a  point  located  a  distance,  b,  radially  out 
from  the  base  of  the  pole.  The  in-line  current  derivative 
sensor  was  located  a  distance,  d,  along  this  line  from  the 
ground  plane  intersection  point.  The  local  angle  of  incidence 
of  the  simulated  EMP  E-field  with  respect  to  the  sensor  is 
indicated  by  a  . 

The  vertical  component  of  the  incident  electric  field 
in  the  vicinity  of  the  current  sensor  on  the  line  is  shown 
in  Figure  3.  This  was  determined  by  measuring  the  horizontal 
component  of  the  magnetic  flux  density,  B  (Figure  3a) ,  and 
then  multiplying  by  the  speed  of  light,  c,  to  obtain  the 
electric  field  (Figure  3b) .  Measurements  of  the  other  field 
components  indicated  that  they  were  at  least  a  factor  of  10 
lower  than  the  vertical  electric  field. 

The  baseline  configuration  chosen  for  comparison  with 
corona  model  predictions  in  this  report  consisted  of  a  86.3 
meter  section  of  1/8"  diameter  aircraft  cable  (7x19  stranding) 


10 


surrounded  by  s  1/2"  PVC  pipe  covered  with  copper  tape 
(OD  =  0.840").  This  was  attached  on  one  end  (A)  to  an 
8  metrr  dielectric  standoff  cable  that  was  then  attached 
to  VPD-II’s  support  pole  #3  (B)  .  The  other  end  was  attached, 
using  standard  fittings,  to  the  OMM-2  current  derivative 
sensor  (C) .  The  output  from  the  sensor  consisted  of  two 
type  141  semi-rigid  coaxial  cables  providing  a  differential 
output,  and  an  aircraft  cable  Lo  provide  a  means  for  tension¬ 
ing  the  entire  test  line  without  stressing  the  coax  output 
cables.  Tension  on  the  test  line  was  achieved  using  a  block 
and  tackle  (D)  attached  to  a  support  "sled"  and  was  monitored 
with  a  dynamometer.  In  general,  the  tension  applied  was  that 
required  to  make  the  line  taut.  Baseline  values  for  the 
variables  indicated  in  Figure  2  were: 

H  =  30.1  m;  L  =  106.3  m;  d  =  12  m;  a  =  18°  . 

As  discussed  in  reference  [2] ,  the  test  line  was  comprised 
of  two  sections.  Above  the  probe,  the  line  was  a  1/8"  diameter 
aircraft  cable.  Below  the  probe,  there  was  a  similar  1/8" 
cable  to  provide  tension  and  two  141  semi-rigid  coaxial  cables 
to  send  the  differential  output  from  the  probe  down  to  a 
fiber-optic  transmitter.  These  three  conductors  ,;ere  then 
covered  with  a  flexible  plastic  tube  and  a  braided  cylindrical 
conductor  was  then  placed  over  the  plastic  tube.  This  resulted 
in  an  equivalent  radius  of  approximately  1.12  cm  for  the 


13 


lower  portion  of  the  test  line.  This  value  is  only  approximate, 
due  to  the  fact  that  the  braided  conductor  did  not  form  a 
perfect  cylindrical  slope  as  it  fit  over  the  plastic  tube. 

In  an  attempt  to  minimise  the  effects  of  having  a  dis¬ 
continuity  in  tha  test  line  radius  at  the  prcbe  location,  two 
different  sizes  of  PVC  tubing  were  covered  with  conducting 
copper  tape  and  placed  over  the  1/8"  aircraft  cable  on  the 
upper  portion  of  the  line.  For  the  nominal  1/2"  PVC  pipe, 
this  increased  the  upper  test  line  diameter  to  0.84"  and  for 
the  1"  PVC  pipe,  the  diameter  was  1.315". 

For  the*  baseline  case  reported  in  [2]  and  analyzed  in 
the  present  report,  the  1/2"  PVC  pipe  was  used  in  an  attempt 
to  minimize  the  discontinuity  of  the  two  radii  of  the  line. 
Figure  4  contains  plots  of  the  data  collected,  both  in 
derivative  and  in  numerically  integrated  forms  fer  the 
baseline  case. 

The  differential  output  from  the  sensor  was  connected 
to  a  Nanofast  300-2A  fiber  optic  transmitter.  Fiber  optic 
cables  transmitted  the  signals  to  the  data  acquisition  system 
(DAS) ,  housed  in  a  shielded  trailer  at  the  far  end  of  the 
VPD-II  working  volume,  where  the  signals  were  recorded  on 
Tektronics  7912AD  waveform  digitizers. 
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Figure  4.  Data  Collected  for  Baseline  Test  Case  During 
Corona  Experiment. 


Tests  involving  other  configurations  similar  to  that 
of  the  baseline  configuration  explored  the  effects  of 
increasing  and  decreasing  the  effective  wire  radius  and 
changes  in  angle  of  E-field  incidence.  Variations  with 
radial  angle  and  tension  were  also  examined.  The  case 
chosen  as  a  baseline  test  for  the  purpose  of  the  comparisons 
in  this  paper  was  selected  for  three  reasons:  (1)  variations 
tested  during  this  experiment  examined  both  increased  and 
decreased  effective  test  wire  radius  as  well  as  angles  of 
incidence  higher  and  lower  than  that  of  the  baseline,  (2)  in 
the  baseline  case,  the  test  wire  radius  was  the  same  as  the 
radius  of  the  OMM-2  probe,  eliminating  concerns  related  to 
possible  discontinuities  at  the  probe  boundary,  and  (3)  the 
diameter  is  representative  of  a  realistic  power  system  con¬ 
ductor.  In  fact,  no  reflections  were  noted  at  that  boundary 
during  any  of  the  testing.  Further  details  of  configuration 
variations  are  contained  in  Reference  [2]  and  [3] . 

In  examining  the  measured  data,  it  is  important  to 
realize  that  there  are  several  perturbations  to  the  response 
which  arrive  at  the  probe  location  at  certain  times.  Figure  5 
shows  the  simulator  and  wire  geometry.  The  first  response 
at  the  probe  is  due  to  a  direct  interaction  of  the  line  with 
the  incident  field  traveling  from  the  pulser  to  the  probe 
along  Path  #1.  Later,  first  order  diffraction  effects  due 
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Direct  Interaction:  Path  #1 

First  Order  Diffraction:  Path  #2 

Path  #3 

Wire  Second  Order  Diffraction:  Path  #4 


Figure  5.  Clear  Times  for  Various  Diffraction 
Components  for  the  Corona  Test. 


17 


to  a  wave  propagating  up  the  simulator  cone  and  then  down 
the  wire  (Path  #2)  and  along  the  ground  and  up  the  wire 
(Path  #3)  will  be  noticed.  Calculations  were  made  prior  to 
the  experiment  to  try  to  maximize  the  clear  time  (i.e.,  the 
time  during  which  the  line  response  is  identical  to  that  of 
an  infinite  line)  and  the  optimum  probe  location  was  determined. 

It  should  be  noted  in  passing  that  there  is  another 
diffraction  component  which  affects  the  response  of  the  line 
current,  although  it  is  of  second-order  importance.  This 
is  shown  as  Path  #4,  which  involves  the  incident  field 
diffracting  from  the  wire,  reflecting  specularly  in  the  earth, 
and  then  impinging  on  the  probe.  This  particular  diffraction 
component  arrives  at  a  time  earlier  than  the  components  from 
Paths  #2  and  #3,  but  is  believed  to  have  a  negligible  effect 
on  the  measured  current.  A  future  research  effort  should 
include  the  determination  of  the  importance  of  this  diffraction 
term. 


During  the  test,  it  was  found  that  the  major  clear  time 
limiting  factor  was  the  reflection  noted  from  the  metallic 
sled  used  in  tensioning  the  line  from  the  ground  end.  In 
general,  the  12  meter  output  from  the  sensor  allowed  sufficient 
clear  time  to  demonstrate  corona  dependent  effects.  Clear 
time  limitations  in  the  measured  data  are  indicated  on 
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each  of  the  data  plots  in  Figure  4,  as  well  as  in  later  plots 
comparing  the  measured  data  with  predictions  made  by  each  of 
the  corona  models  studied  in  this  report. 

B.  Data  Processing  Concepts 

The  differential  output  from  the  OMM-2  probe  was  attached 
to  a  fiber  optic  transmitter  containing  a  balun  for  common 
mode  rejection.  The  measured  signal  was  then  sent  over  fiber 
optic  cables  to  the  data  acquisition  system,  housed  nearby  in 
a  shielded  trailer.  At  that  point,  the  signal  was  temporally 
split,  using  an  active  signal  divider,  and  passed  to  waveform 
digitizers  such  that  there  was  a  5  to  10  nanosecond 
overlap  of  data  in  successive  digitizers.  This  overlap 
allowed  the  data  to  be  later  aligned  and  connected  on  a  common 
time  scale  (a  process  called  "time-tying"  the  data) .  After 
the  data  was  recorded  on  the  waveform  digitizers,  it  was 
transferred  to  the  TRESTLE  VAX  /  11-750  computer  for 
subsequent  processing  by  BDM. 

The  data  processing  techniques  discussed  below  apply  to  a 
collection  of  data  sets  representing  the  line's  response  to  a 
single  VPD  shot.  The  term  "subsequent  data  set"  means  the  set 
of  data  recorded  by  a  single  waveform  digitizer  that  is 


19 


temporally  next  to  (and  overlapping  with)  the  data  set  just 
considered. 

On  the  VAX  11-750  ,  the  data  was  linearly  scaled  to 
correct  for  non-frequency  dependent  instrumentation  effects. 
Data  from  the  first  digitizer  was  also  corrected  for  any 
obvious  DC  offset.  The  assumption  made  in  the  DC  offset 
correction  is  that  in  the  absence  of  any  excitation  the  mean 
of  the  line  response  data  should  be  zero. 

Next,  data  from  subsequent  digitizers  is  time-tied  wi\.i- 
that  from  earlier  digitizers.  Beginning  with  the  tail  of  the 
first  set  of  data  and  the  head  of  the  second  set  of  data, 
contemporaneous  sections  of  data  are  compared  and  blended 
through  averaging  and  offsetting  the  data  as  necessary.  Note 
that  the  DC  correction,  applied  only  to  the  earliest  data  set, 
is  passed  through  subsequent  data  sets  without  requiring  a 
priori  knowledge  or  correction  of  the  exact  DC  offset  of  those 
later  data  sets.  This  DC  offset  adjusted,  scaled,  and 
time-tied  data  is  referred  to  as  the  raw  data. 

Raw  data  may  still  contain  some  frequency  dependent 
instrumentation  bias.  If  the  transfer  function  of  the 
combined  effects  of  all  instrumentation  does  not  differ 
significantly  from  a  constant  value,  its  deconvolution  may  be 
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performed  in  the  time  domain  by  further  scaling  of  the  data  by 
the  amount  of  that  constant.  However,  if  the  frequency 
dependencies  are  critical,  it  is  easiest  to  deconvolve  them  in 
the  frequency  domain.  Here,  a  Fourier  integral  transform  was 
used  to  convert  the  time  domain  data  to  the  frequency  domain 
and  back  to  the  time  domain  following  deconvolution  of 
instrument  effects.  Frequency  domain  transfer  functions  for 
each  indicated  piece  of  equipment  consisted  of  50  points  per 
decade  over  the  three  decade  range  250  kHz  to  250  MHz  and 
were  obtained  using  a  Hewlett  Packard  8505A  Network  Analyzer 
calibrated  to  NBS  standards  by  Precision  Measurement  Equipment 
Laboratory  (PMEL) . 

C.  Measurement  Equipment 

Table  I  contains  a  list  of  the  major  equipment  used  to 
conduct  this  experiment. 


TEST  LINE 

1/8”  (diameter)  aircraft  cable  (7x19  stranding) . 
1/8”  aircraft  cable  with  1/32”  thick  PVC  coating. 
1/2"  (inner  diameter)  schedule  315  PVC  pipe. 

1"  (inner  diameter)  schedule  200  PVC  pipe. 

Copper  tape  added  to  exterior  of  all  PVC  pipes. 
EG&G  OMM-2  Differential  current  probe 
EG&G  DMB-1  (differential  mode  balun) 

Tinned  copper  grounding  braid. 

DATA  ACQUISITION 

Tektronics  7912AD  Waveform  Digitizers  with  plug  in 
modules : 

Tek  7A16P  Programmable  Amplifier 
Tek  7B90P  Programmable  Time  Base 
Nanofast  300-2A  Fiber  Optic  Link  -  Transmitter 
Hewlett-Packard  11549A  Power  Splitter 
Tektronics  7104  1  GHz  Real  Time  Oscilloscope  & 
modules : 

Tek  7A29  Dual  Trace  Amplifier 
Tek  7B15  Delta  Delaying  Time  Base 
EG&G  ODL-6  Fiber  Optic  Link  -  Receiver 

MISCELLANEOUS 

Omni-Spectra  2082-6013-20  20  dB  attenuators 
EG&G  MGL-2  B-dot  sensor  (for  field  mapping) 


Table  I. 


Major  Equipment  used  in  Corona  Experiment 


D.  Limitations  of  the  Experiment 


It  should  be  noted  at  this  point  that  there  were  limits 
to  what  could  be  examined  in  the  course  of  the  experiment. 
Furthermore,  the  analysis  contained  below  is  limited  since 
it  examines  the  predictions  of  each  analytic  model  against  a 
single  baseline  test  case,  selected  from  the  full  set  of 
acquired  data  for  the  reasons  discussed  earlier. 

Consequently,  the  reader  may  be  aware  of  two  kinds  of 
limitations:  those  arising  from  limits  in  what  could  be 
tested  during  the  course  of  the  experiment  and  those  arising 
from  the  use  of  a  single,  representative  baseline  test  case  in 
evaluating  the  models  discussed  below.  Such  limitations  are 
inherent  in  nearly  every  experiment  and  in  nearly  every 
comparison  of  analytical  results  with  experiment  data. 

Limitations  on  the  experiment  included  limitations 
posed  by  the  capabilities  of  the  EMP  simulator  used,  as  well 
as  those  created  by  the  geometry  chosen  ft>r  the  experiment. 

The  Incident  field  effectively  had  a  fixed  form  with  a 
constant  rise  time.  Consequently,  the  induced  current  had 
essentially  a  constant  rate  of  rise.  Different  waveforms  may 
have  an  impact  on  the  corona  incepcion  voltage,  but  there  is 
no  data  from  this  experiment  that  clearly  supports  or  negates 
this  notion.  Given  the  full  set  of  acquired  data  [2,  3],  one 
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night  be  able  to  make  tentative  conclusions  about  the  corona 
inception  voltage,  but  such  conclusions  would  be  inappropriate 
in  this  paper  due  to  the  United  nature  of  the  single  baseline 
test  case.  However,  sinulator  linitations  restricted  the 
variations  in  the  pulse  shape  used  during  the  experinent. 
Consequently,  even  the  full  set  of  acquired  data  nay  be 
inadequate  for  full  detemination  of  the  relation  of  the 
corona  inception  voltage  to  the  pulse  risetine. 

There  are  two  polarities  for  the  EMP.that  could  be 
tested  to  determine  the  difference  in  their  inpact.  However, 
the  pulser  of  the  sinulator  could  only  be  charged  in  one  way, 
resulting  in  an  E-field  that  generated  a  positively  charged 

wire.  In  addition,  energizing  the  test  line  (to  give  it 

\ 

electrical  characteristics  similar  to  either  an  AC  or  a  DC 
transmission  line)  night  also  impact  the  results  noted. 

Again,  we  have  no  data  to  support  or  negate  these  ideas.  They 
are  presented  here  to  emphasize  the  linitations  of  what  might 
be  inferred  from  the  data. 

Reflections  of  an  EMP  from  the  ground  plane  such  that 
its  initial  interaction  with  the  test  line  is  followed  60-100 
nanoseconds  later  by  the  interaction  of  the  reflected  EMP  with 
the  test  line  were  not  tested  because  both  the  simulator 
and  the  experiment  geometry  could  not  accomodate  such  a  test. 
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The  reflected  pulse,  which  is  likely  to  occur  as  a  result  of 
an  EMP  resulting  from  a  nuclear  weapon  detonation,  can  be 
expected  to  also  alter  the  results  obtained.  However  it  will 
be  noted  that  none  of  the  models  investigated  here  considered 
this  reflected  EMP  either.  Although  the  reflection  is  likely 
to  occur  in  an  actual  EMP  event,  the  magnitude  of  its  impact 
on  a  given  transmission  line  is  strictly  a  function  of  the 
local  geometry.  Later  revisions  to  corona  models  n iy  include 
considerations  of  the  impact  of  the  pulse  reflection  in  terms 
of  the  local  geometry.  The  consideration  of  the  reflected  EMP 
is  beyond  the  scope  of  the  current  investigation. 

Although  a  number  of  attempts  were  made  during  the 
experiment  to  obtain  "no  corona"  data  [2,3],  the  effect  of 
corona  was  noted  in  every  test  run.  However,  in  some  of  these 
attempts,  the  reduction  in  corona  impact  was  sufficient  to 
indicate  a  general  notion  of  what  "no  corona"  data  should  look 
like.  A  comparison  of  the  baseline  test  data  with  analytical 
"no  corona"  data  is  included  in  Figures  12-20. 
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IV.  COMPARISON  OF  MEASURED  DATA  WITH  CALCULATIONAL  RESULTS 

A.  Overview 

The  study  of  the  effects  of  corona  discharge  on  power 
system  conductors  has  been  a  subject  of  considerable  interest 
for  a  number  of  years.  This  interest  stems  from  corona's  role 
in  creating  radio  frequency  interference,  in  adding  power  loss 
on  electrical  transmission  lines,  and  in  modifying  the  shape 
of  lightning- induced  transient  surges  along  power  lines.  A 
number  of  authors  [5],  [10],  [11],  and  [12]  have  studied 
corona  related  effects  which  are  frequently  described  as 
arising  from  a  changing  conductor  capacitance  ar  ’  a  corona 
resistance  which  accounts  for  loss  in  the  corona. 

In  the  recent  literature,  there  are  a  number  of  articles 
which  describe  corona  models  for  power  lines  or  other 
above-ground  conductors  [13],  [14],  and  [15].  All  of  these 
models  make  the  assumption  that  the  propagation  of  voltage  and 
current  waves  on  the  line  can  be  described  by  TEM  transmission 
line  theory,  with  suitable  modifications  being  made  in  either 
the  source  terms  or  in  the  transmission  line  parameters  to 
account  for  the  behavior  of  the  corona. 
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Although  this  previous  work  aided  in  the  development 
of  corona  models,  its  application  was  primarily  intended  for 
the  study  of  point-injected  signals,  such  as  that  form  a 
lightning  attachment  to  the  line,  and  their  subsequent  propa¬ 
gation  characteristics.  Other  early  corona  models  examined 
corona  resulting  from  high  voltages  being  carried  on  the  lines. 
Neither  of  these  forms  of  corona  model  were  directly  applicable 
to  the  case  of  EMP  induced  corona. 

In  the  EMP  area,  it  is  important  to  understand  the 
behavior  of  a  distributed  field  excitation  of  a  line  as  well 
as  the  effects  that  EMP-induced  corona  will  have  on  the  line 
current.  There  have  been  a  number  of  recent  papers  which 
describe  models  for  computing  this  phenomena  [6],  [7],  [16], 

[17]  and  [18],  As  in  the  previous  cases  of  surge  propagation 
in  the  presence  of  corona,  most  of  these  models  invoke  the 
transmission  .line  approximation  for  obtaining  a  solution  for 
the  line  response,  but  with  a  distributed  field  excitation. 
Investigations  have  indicated  that  such  transmission  line 
solutions  can  be  quite  accurate  in  describing  the  transient 
current  behavior  as  illustrated  in  Figure  6. 

An  alternate,  but  closely  related  calculational  approach 
is  described  in  reference  [7]  which  employs  an  antenna  solution 
for  the  corona  coupling  problem.  This  and  che  other  trans¬ 
mission  line  coupling  models  will  be  described  in  more  detail 
in  the  following  section. 
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ct/th  an 9) 


Figure  6.  Current  Induced  on  an  Infinite  Conductor 
Over  a  Perfectly  Conducting  Ground  Plane 
by  a  Step-Function  Incident  Wave  (Trans¬ 
verse  Magnetic  Field)  [25] . 


In  the  corona  models  to  be  discussed,  there  are  several 
parameters  which  describe  the  corona  behavior  and  its  effect 
on  the  line  response.  Frequently,  the  precise  values  for  these 
corona  model  parameters  are  not  known  from  first  principles. 
However,  it  is  believed  that  they  can  be  determined  by  com¬ 
paring  computed  data  with  experimental  results. 

One  common  parameter  in  these  models  is  the  air  breakdown 

electric  field  strength,  Ec,  which  is  usually  taken  to  be  a 

constant.  Ec  is  considered  to  be  independent  of  the  rate  of 

rise  of  the  local  electric  field  normal  to  the  conductor, 

which  is  the  ultimate  ionization  field  for  the  air  surrounding 

the  conductor.  It  may  be  noted  that  in  several  of  the  models 

studied,  better  fits  to  the  data  have  been  found  using  values 

for  Ec  larger  than  the  range  of  values  predicted  by  Peek's 

law.  Possible  reasons  for  this  may  include:  (1)  There  is 

a  lag  between  the  time  when  the  local  field  equals  or  exceeds 

E  and  the  time  when  the  impact  of  corona  formation  is  noted 
c 

in  the  conductor  response.  This  time  is  required  for  appropriate 
numbers  of  separated  electron-ion  pairs  to  form.  If  the 
driving  waveform  has  a  magnitude  higher  than  E  and  risetime 
slightly  longer  than  the  corona  establishment  time,  then  the 
time  delay  in  corona  establishment  could  be  seen  in  the  data 
as  an.  artificially  high  Ec>  Reference  4  indicates  that  the 
establishment  time  consists  of  a  statistical  lag  time,  which 
is  dependent  on  geometry  and  gas  characteristics,  and  a 
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formative  lag  time  (about  1  nanosecond) .  In  this  particular 
experiment,  the  statistical  lag  time  depends  on  the  time  for 
the  electron  to  detach  from  an  O2  ion  in  the  region  in  which 
the  electric  field  exceeds  E  [20],  (2)  Peek's  law  as  stated 

v 

in  [4]  is  valid  only  from  DC  to  several  kiohertz.  The  driving 
waveform  contains  higher  frequencies,  although  the  bulk  of  its 
energy  is  in  the  lower  frequencies.  The  impact  of  the  higher 
frequency  components  may  be  to  raise  the  effective  E  ,  although 
there  is  evidence  (19,  p.  4-136]  contrary  to  this  notion. 

In  this  section,  the  measured  EMP-induced  response  of  a 
wire  excited  strongly  enough  to  produce  a  corona  along  the 
conductor  is  compared  with  5  different  corona  models.  All  of 
there  models  involve  corona  parameters  which  must  be  "fit"  to 
a  particular  set  of  measured  data,  then  used  to  predict  the 
response  of  the  line  to  other  excitations. 

This  section  is  divided  into  several  sub-sections.  Each 
of  these  will  describe  one  of  the  previously  referenced  corona 
models  and  compare  the  corresponding  computed  line  responses 
with  the  measured  results  which  are  described  in  refernece  [2] . 
From  this  study,  a  suggestion  as  to  the  "best"  corona  model 
for  BMP  predicitons  can  be  made. 
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B.  The  Baum  Corona  Model 


1 .  Theory 


The  model  proposed  by  Baum  [6]  for  the  EMP  excitation  of 
a  conductor  is  similar  to  previous  corona  models  used  within 
the  power  community  for  computing  the  propagation  of  lightning 
surges  along  a  line.  However,  in  the  EMP  case,  a  set  of 
distributed  excitation  sources  along  the  line  is  required.  In 
this  model,  a  voltage-dependent  conductor  radius  is  used  to 
represent  the  effects  of  the  corona,  and  the  conductor  is 
assumed  to  be  located  in  free  space  or  above  a  reference 
ground  plane  at  a  height,  h  (Note  that  Baum's  model  describes 
a  conversion  between  a  conducting  line  above  a  ground  plane 
and  an  analogous  coaxial  line.).  The  response  of  the  line  is 
computed  by  transmission  line  equations  in  which  the  line 
current  and  voltage  (with  respect  to  the  ground  plane)  are 
unknowns : 


av 

dz 


+ 


V' 


s 


(la) 


and 
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(lb) 


91 

dz 


+ 


|t(C-V) 


In  this  expression,  I(z,t)  is  the  induced  current  on 

the  line  and  V(z,t)  is  the  voltage  of  the  line.  The  term  L' 

is  the  per-unit-length  inductance  of  the  line,  and  this  is 

assumed  to  be  independent  of  the  presence  of  corona  on  the 

line.  For  the  case  of  a  line  of  radius  a„  located  over  a 

w 

ground  at  a  height  h  ,  this  inductance  is  given  by 


L‘  =  ^o  fgo  ' 


where  the  parameter  f  is  given  by 


f  .  =  1  cosh-1  (h/arr) 

9°  2tT~  w 


2t r 


In  (2h/aw) 


(2) 

(3a) 

(3b) 


The  term  C*  is  the  per  unit  length  capacitance  of  the 
line.  In  the  Baum  model,  this  is  assumed  to  depend  on  the 
nature  of  the  corona.  Normally,  this  parameter  is  given  by 


C* 


Vfgo 


(4) 


for  the  case  of  a  wire  without  corona.  With  corona,  however. 
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Baum's  model  assumes  that  the  wire  radius  is  effectively 
increased,  due  to  the  presence  of  the  ionized  cloud 
surrounding  the  wire.  Consequently,  the  per  unit  length 
capacitance  is  modified.  The  new  capacitance  is  given  by 


C' 


eo/fgc 


(5) 


where  f  has  the  same  form  as  in  equation  (3) ,  with  the 
radius  of  the  wire,  a  ,  being  replaced  by  the  effective 

W 

radius  of  the  corona  region,  ac  . 

This  corona  radius  is  estimated  by  determining  the 
radius  away  from  the  wire  at  which  the  local  electric  field 
just  exceeds  the  air  breakdown  field  strength.  Denoting  this 
breakdown  field  by  E  ,  the  radius  of  the  corona  region  can 
thus  be  expressed  in  terms  of  the  charge  per-unit-length  on 
the  conducting  wire,  Q'  ,  as 


2  v  &  E_ 
o  c 

The  line  capacitance  is  given,  therefore,  by  the 
nonlinear  relation 
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(7) 


C'(Q)  =  £0/ (minimum  [fgo, fgc(Q) ] )  , 

which  depends  on  the  local  charge  on  the  line,  as  well  as  on 
the  air  breakdown  strength. 

For  the  case  of  an  isolated  wire  in  free  space,  the  i-v 
relations  of  equations  (la)  and  (lb)  are  still  assumed  to 
hold,  but  the  definition  of  the  line  inductance  in  equation 
(2)  and  the  line  capacitance  in  equation  (4)  becomes 
ill-defined,  due  to  the  lack  of  a  reference  conductor  in  the 
problem.  In  this  instance,  the  parameter  f  (or 
equivalently,  the  line  "height”,  h)  may  be  regarded  as  a 
fictitious  line  constant  and  varied  in  order  to  obtain  a  good 
fit  of  the  computed  data  with  that  obtained  from  experiment. 
Note  that  the  response  of  the  line  depends  on  h  only 
logarithmically,  so  as  a  consequence,  the  response  is  not 
extremely  sensitive  to  the  choice  of  this  parameter. 

The  source  terms  in  equations  (la)  and  (lb)  are  related 

to  the  incident  field  exciting  the  line.  The  term  is 

related  to  the  per-unit-length  magnetic  flux  linking  the 

conductor  and  the  reference  conductor,  and  the  term  I'  is 

s 

related  to  the  charge  deposited  on  the  conductor  due  to  the 
incident  electric  field. 
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In  reference  [6],  Baum  illustrates  how  these  equations 
can  be  re-arranged  to  yield  two  coupled  equations  for  the  line 
current  and  charge,  with  the  source  term  related  to  the 
incident  tangential  electric  field  on  the  wire,  E^.  .  This 
form  is  completely  equivalent  to  the  previous  equations,  and 
is  expressed  as 


fl  (Q'/C)  +  L' 


(8a) 


and 


31  9Q*  _ 

dz  at  “ 


0 


(8b) 


By  defining  a  new  variable  v(z,t)  s  Q(v,t)/c’ (Q)  ,  it 
is  possible  to  convert  equations  (8)  into  the  following 
"standard"  form,  which  will  be  useful  in  comparing  this  corona 
modeling  approach  with  other  models: 


dv 

dz 


L.  31  _ 

L  at  " 


and 


(9a) 


31 

dz 


+ 


§t<c'v)  -  0 


(9b) 
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Note  that  although  the  variable  v  has  the  dimensions 
of  volts,  it  is  not  the  total  voltage  between  the  line  and  the 
ground  which  is  being  induced  by  the  EMP.  This  latter  response 
is  computed  by  adding  the  contribution  of  the  vertical 
component  of  the  incident  electric  field  as  described  in  [4]. 

It  is  the  set  of  equations  in  (9)  that  must  ultimately 
be  solved  for  the  line  current,  I  .  This  solution  may  be 
carried  out  using  the  usual  time-marching,  finite-difference 
approach,  or  by  a  method  reminiscent  of  the  method  of 
characteristics  [21],  as  done  by  Baum. 

In  applying  this  corona  model  to  the  experimental  data 

described  in  Section  II,  it  is  noted  that  there  are  two 

independent  parameters  which  can  be  varied  in  order  to 

optimize  the  computed  solution.  These  are  the  electric  field 

breakdown  strength  in  air  E„  ,  and  the  value  of  f _  in 

3  c  go 

equation  (3)  which  can  be  interpreted  as  fixing  the  value  of 
the  parameter  h  . 

Considering  the  geometry  for  the  baseline  experimental 

case  as  defined  in  Section  II,  corona  coupling  calculations 

using  the  Baum  model  have  been  conducted  for  a  variety  of 

different  values  of  E„  and  f _  ,  and  the  results  are 

c  go 

presented  in  figures  7  through  20. 
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2.  Discussion  of  Results 


In  this  model  there  are  several  parameters  which  are  not 
well  known  from  the  geometry  of  the  test  configuration,  and 
hence,  may  be  varied  in  order  to  obtain  a  reasonable  fit  of 
the  computed  data  to  the  measured  results.  These  variables  and 
the  computed  results  are  discussed  in  more  detail  below. 

This  subsection  is  comprised  of  the  following  elements: 
a  list  of  the  model  variables  with  a  brief  description  of  each 
variable;  known  limitations  on  the  range  of  each  variable;  a 
list  of  the  specific  values  of  each  variable  tested;  and  a 
results  summary  in  which  the  results  are  described  in  terms  of 
how  changes  in  each  variable  affect  the  fit  of  the  model  to 
the  data  and  how  well,  in  general,  the  model  could  predict  the 
results  acquired  experimentally. 

Variables: 


E  E-field  value  at  which  breakdown  occurs 
c 

f  dimensionless  factor  to  convert  line  geometry  to 

equivalent  coaxial  transmission  line  or  line  above 
conducting  plane 
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Limitations : 


6  6 

Peek's  law  constrains  2.5  x  10  V/m  <  Ec  <  4.5  x  10  V/m. 

As  indicated  earlier,  for  a  line  in  free  space  above  a 

conducting  ground  plane,  f  =  1/2tt  *  In  (2h/aw)  .  For  a 

coaxial  transmission  line,  f _ =  l/2ir  *  In  (r  .  ^ /r.  . 

go  outer'  inner 

In  the  experiment,  the  test  wire  ran  diagonally  from 
approximately  30  meters  above  ground  to  about  1  meter 
above  the  ground  plane.  This  geometry  does  not  readily 
translate  into  a  constant  value  for  f  ,  although  for  this 

range  of  wire  height,  the  range  of  f  is  fairly  small: 

1.14  to  1.68.  In  general,  realistic  values  of  f  lie  in 
the  range  0.1  <  f go  <5.0  . 

Variations  tested: 

E  :  2.5  X  106,  3.0  X  106,  3.5  X  106,  4.0  x  106, 

c 

4.5  x  106  (all  in  V/m) 

f  :  0.5,  0.75,  1.0,  1.2,  1.4,  1.6,  1.8,  2.0,  3.0 
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Results: 


Figures  7-11  show  the  effect  of  changing  f  while 

holding  Ec  constant.  Essentially,  higher  values  for  f 

reduce  the  initial  slope  of  the  curve  (i.e.  the  "no  corona" 

section)  and  reduce  the  magnitude  of  the  corona  suppressed 

current.  In  all  5  figures,  the  curve  for  f  =  1.2 

appears  to  be  the  best  fit.  This  corresponds  to  an  equivalent 

geometry  in  which  the  test  line  would  be  approximately  10 

meters  above  the  ground  plane.  Alternatively,  this 

corresponds  to  a  equivalent  coaxial  cable  with  a  ratio  of 

r  .  to  r. _  of  about  1880  ,  or  a  coaxial  cable  with 

outer  inner 

the  test  line  used  as  the  center  conductor  (r  =  1.07  cm) 
surrounded  by  an  outer  conductor  of  radius  20  meters. 

Figures  12-20  indicate  the  effect  of  changing  E  while 

holding  fgQ  constant.  Note  that  increasing  values  for  Ec 

demonstrate  delayed  corona  onset  times,  as  expected.  Figure 

15  (showing  variations  in  E  for  a  constant  f  =1.2 

c  go 

6 

indicates  E  =  5.25  x  10  V/m  is  approximately  the  value  at 
c 

which  this  model's  prediction  of  corona  onset  matches  the 
experimental  data.  It  is  possible  that  this  value  may  be 
artificially  high,  since  corona  onset  is  not  instantaneous,  as 
this  model  assumes.  In  fact,  corona  is  a  complex  function  of 
air  chemistry.  It  may  be  initiated  earlier  in  time  at  a  lower 
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nummiwue 


E  ,  but  its  full  impact  on  the  suppression  of  the  line 
c 

current  may  not  be  seen  until  the  corona  is  fully  developed 
(i.e.,  until  a  sufficient  amount  of  the  air  surrounding  the 
test  line  is  ionized,  essentially  forming  a  capacitance 
between  the  line  and  the  corona  sheath) . 

In  general,  the  model  fit  to  the  experimental  data  is 
fair.  The  model  prediction  of  current  suppression  is  less 
than  actually  noted.  This  may  be  due  to  several  factors:  (1) 
the  model  assumption  of  instantaneous  corona  formation 
neglects  the  current  suppression  effects  during  corona 
formation,  (2)  the  capacitance  assumed  between  the  corona 
sheath  and  the  test  line  may  be  too  small,  (3)  the  model 
assumes  a  perfectly  conducting  test  line.  A  revision  of  the 
model  to  include  inherent  line  losses  might  bring  the 
predictions  more  in  line  with  those  seen  experimentally. 
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C.  The  Engheta-Lee  Corona  Models 


As  previously  mentioned,  reference  [7]  takes  a  slightly 
different  approach  for  the  modeling  of  a  line  with  corona. 
Starting  from  Maxwell's  equations  and  the  boundary  conditions 
on  the  wire,  the  following  expressions  for  current  and  charge 
on  the  surface  of  the  conductor  are  derived: 


1  3£'  31  _  aE1 

~  ”  +  n  —  =  2ira.  -5— 

e  dz  o  3t  w  dp 

o 


(10a) 


and 


dz  at 


-  2ira  Jn 
w  p 


(10b) 


These  equations  are  identical  to  those  developed  in  [15], 


With  the  definition  that  v(z,t)  =  Q’ (z,t)/eo  ,  this 
equation  can  be  put  into  the  standard  form  given  in  equation 
(9)  as: 


av 

dz 


31  ,  5Et 

^O  at  =  2iraw  Wp 


(Ha) 


and 
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(lib) 


51 

5z 


+ 


a_ 

at 


(e0v) 


2ira  J 

P 


Note  that  in  this  formulation  there  is  no  ground  plane 

introduced.  In  comparing  equations  (11)  with  the  previous 

equations  (9)  it  is  seen  that  the  per-unit-length  inductance 

and  capacitance  of  the  line  have  been  replaced  by  the 

permeability  and  permittivity  of  free  space,  and  that  there  is 

a  source  term  in  equation  (lib).  This  source  is  proportional 

to  J  and  in  the  Engheta-Lee  paper,  two  models  are  discussed 
P 

in  which  Jp  takes  one  of  two  possible  forms:  one  arising 
from  considering  the  corona  as  a  lossy  conductor  (conductivity 
model),  and  the  other  from  considerations  of  Townsend's  theory 
of  corona  (Townsend  model)  [8].  Both  of  these  source  terms 
will  be  discussed  in  more  detail  in  this  sub-section. 


The  other  source  term  in  equation  (11a)  is  related  to 
the  radial  derivative  of  the  total  electric  field  in  the  z 
direction.  Unfortunately,  the  behavior  of  this  component  of 
the  field  is  unknown,  and  not  easily  related  to  the  incident 
electric  field.  In  fact,  this  difficulty  was  discussed  in  [15] 
and  resulted  in  a  modified  approach  for  the  solution  of  the 
line  response  in  that  reference. 

In  the  development  in  [7],  this  troublesome  term  is 
determined  by  examining  the  behavior  of  a  single  conductor  in 
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free  space  without  corona.  For  this  case,  reference  [22] 

develops  a  set  of  equations  similar  to  those  of  (11),  and  it 

is  noted  that  the  current  source  term,  J  ,  is  zero  and  the 

P 

source  term  for  the  first  equation  is  related  to  the 
tangential  component  of  the  incident  electric  field.  In  the 
Engheta-Lee  corona  models,  it  is  assumed  that  the  excitation 
term  in  equation  (11a)  can  be  replaced  with  that  arising  from 
the  incident  field,  and  that  the  two  telegrapher's  equations 
for  the  line  response  then  take  the  following  form: 


av  ai 

3z  ^o  3t 


4w 

n(t) 


E“°(t) 


(12a) 


and 


~  +  I ~(e  v)  =  -  2rra„  J  ,  (12b) 

dz  3tv  O  '  W  p  '  v  ' 

where  the  factor  fi  is  a  time-varying  quantity  given  by 

Q(t)  =  2  ln[  2 (ct  -  zcos(0)  +  a  sin(9))  /  (Ta  sin(0))  ]  (13) 

W  W 

with  T  =  exponential  of  Euler's  constant  =  1.7810...  and 

0  =  angle  of  E-field  incidence  with  respect  to  the  wire. 
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Notice  that  doing  this  assumes  that  the  total  field 

term,  OE^/dp,  does  not  depend  on  and  is  not  altered  by  the 

presence  of  the  corona  around  the  wire.  The  validity  of  this 

and  other  assumptions  regarding  the  form  of  J  will  be 

P 

ultimately  validated  in  comparing  the  experimental  response  of 
the  line  with  the  data  calculated  as  discussed  below. 

1.  The  Conductivity  Model 

a .  Theory 

In  the  conductivity  model  postulated  by  Engheta  and  Lee, 
it  is  assumed  that  when  the  line  is  driven  into  corona,  the 
radial  component  of  the  corona  current  density  Jp  is 
proportional  to  the  radial  electric  field  Ep  at  the  wire 
surface  a* 

J  =  o  E  U(|E  I  -  E_)  ,  (14) 

p  p  x  1  p  1  c#  # 

where  E„  is  the  critical  air  breakdown  electric  field 
strength,  a  is  the  effective  conductivity  of  the  corona  and 
U(  )  is  the  Heaviside  unit  step  function.  This  step  function 
assures  that  the  conductivity  effects  of  the  corona  are  not 
seen  in  the  model  if  the  local  field  strength  on  the  wire  is 
below  the  critical  air  breakdown  field  strength. 
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This  last  expression  involving  the  electric  field  may  be 
related  to  the  local  charge  density  on  the  wire  and  can  be 
rewritten  as 


,'r  U  [  SaV  -  Ec  ]  • 


In  terms  of  the  variable  v  =  Q‘/e0  »  this  expression  for  the 
source  term  becomes 


2lfa  U  (  2ira  “  Ec  ) 


The  conductivity  model,  therefore,  requires  a  solution 
of  equations  (12a)  and  (12b)  with  Jp  being  determined  by 
equation  (16) .  It  involves  the  specification  of  two 
independent  parameters:  the  critical  electric  field  strength 
Ec  and  the  corona  conductivity  a  .  The  method  of 
characteristics  [21]  was  used  in  [7]  to  obtain  a  solution, 
however,  in  the  present  study,  a  direct  time-marching  solution 
of  equations  (12a)  and  (12b)  was  used.  Comparisons  between  the 
responses  using  the  two  calculational  methods  yielded 
identical  results. 


For  this  model,  a  parametric  study  was  performed  using 

the  baseline  geometry,  in  order  to  ascertain  the  best  choice 

for  E  and  a  .  Figures  21  through  35  present  the  results  of 
c 

these  calculations. 
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21 .  Comparison  of  Measured  ( - )  and  Conductivity 

Model  ( - )  Results: 

Ec  =  1.0  MV/m;  o  Variations  as  Indicated 


a(mhos/m) 


o,  10" 


10 


-5 


10 


-4 


10 


10 


-2 


22.  Comparison  of  Measured  ( - )  and  Conductivity 

Model  ( - )  Results: 


„  -  1.5  MV/m;  o  Variations  as  Indicated 


E 


2.00E+3 


Figure  31 .  Comparison  of  Measured  ( - )  and  Conductivity 

Model  ( - )  Results: 

o  =  5x10  5  mhos/m;  Variations  as  Indicated 

c 


Figure  32.  Comparison  of  Measured  ( - )  and  Conductivity 

Model  ( - )  Results: 

-4 

=  10  mhos/m;  E_  Variations  as  Indicated 


o 


b.  Discussion  of  Results 


The  following  subsection  discusses  the  limits  of  the 
parameters  of  this  corona  model  and  the  resulting  computed 
responses . 

Variables; 

E  E-field  value  at  which  breakdown  occurs 
c 

a  air  conductivity 

Limitations: 

Peek's  law  constrains  2.5  x  106  V/m  <  Ec  <  4.5  x  106  V/m. 

Air  conductivity  during  corona  is  very  difficult  to 
measure  since  it  is  a  function  of  the  amount  and  polarity  of 
ionization  as  well  as  of  the  local  relative  humidity  and 
pressure  (alternatively,  it  is  a  function  of  the  mobility  and 
number  density  of  the  charge  carriers) .  Reference  6  suggests 
some  possible  values  for  use  with  this  model.  For  comparison, 
some  nominal  conductivity  values  for  well  known  materials  are 
given  in  Table  II  below  [taken  from  reference  22]. 


Table  II. 

Nominal  Conductivity 
Values  for  Selected 
Materials 


Material 

Conductivity.  e{ S/m) 

Material 

Conductivity.  n\ S  ml 

Silver 

6.n  *  io’ 

Freeh  we  ter 

10‘ 5 

Copper 

5.80  V  10’ 

Distilled  wetcr 

2  x  IO'4 

Gold 

4.10  x  10’ 

Dry  toil 

10-’ 

Aluminum 

3.54  x  10’ 

Transformer  oil 

I0-" 

Brass 

1.57  x  10’ 

Glue 

to*1 

Bronze 

10’ 

Porcelain 

2  x  ter” 

iron 

10" 

Rubber 

10” 

Seawiter 

J 

[  Fueed  quart! 

IJ1’ 
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Variations  tested: 


Ec:  1.0  x  106,  1.5  X  106,  2.0  x  106,  2.5  x  106, 

3.0  x  106,  3.5  x  106  (all  in  V/m) 

a:  10~6,  10-5,  10~4,  10-3,  10-2  (all  in  mhos/m) 

Results : 

Figures  21-28  show  the  predictions  of  this  model  for 

constant  values  of  E  .  For  larger  values  of  E  ,  the  model 

c  c 

predictions  begin  to  approximate  measured  current  values  when 
—4 

a  -  10  mhos/meter. 

Holding  a  constant  while  varying  E  is  demonstrated 

C 

in  Figures  29-35.  Extremely  low  values  of  a 

(  or  <  l0~6n-1m-1}  show  no  variation  with  E„  .  This  is  to  be 

c 

expected,  since  for  those  values  of  a  ,  the  air  remains  an 

insulator  at  even  high  field  strengths.  For  a  -  10~5n”1m"’1 

some  variation  is  seen  with  Ec  ,  but  the  model  fit  to  the 

experimental  data  is  still  not  especially  good.  At  a  = 

-4  -1  -l 

10  fl  m  ,  the  results  predicted  by  the  model  are  quite 
close  to  those  seen  experimentally,  especially  for  higher 
values  of  E_  .  At  higher  values  of  a  ,  the  model  predictions 
become  irregular,  possibly  due  to  instabilities  in  the  method 
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of  solution  used.  These  are  especially  noticeable  at  a  = 

-3  -1  -1  -3  -1  -1 

10  n  m  .  Values  of  a  above  10  n  m  produce  results 

well  below  the  observed  values. 

2 .  The  Townsend  Model 


a .  Theory 

The  Townsend  model  as  described  in  [7]  assumes  that  the 
current  term  in  equation  (12b)  has  the  form 

J  =  a.  Q‘  f  - E  1  U  f  -■--2  J - E  1  (17) 

P  J  l  2iraveo  c  J  l  2iraweo  c  j 

or,  in  terms  of  the  variable  v, 


J 

P 


=  a  .e 
J  ° 


(18) 


The  term  in  the  above  equations  is  related  to  a 

Townsend  ionization  coefficient  as  discussed  in  [7],  and  in 
this  model  is  a  parameter  to  be  varied  to  fit  the  experimental 
data.  As  in  the  previous  case,  a  series  of  calculations  for 


the  baseline  line  geometry  was  made  with  a.  and  E  varying, 

D  c 

and  these  results  are  presented  in  Figures  36  through  48 . 
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Figure  36.  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 

E  =1.0  MV/m;  a.  Variations  as  Indicated 
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Figure  37.  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 

E  =1.5  MV/m;  a.  Variations  as  Indicated 
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Figure  38.  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 

E  *  2.5  MV/m;  a .  Variations  as  Indicated 
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Figure  39.  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 

E  =2.0  MV/m;  a.  Variations  as  Indicated 
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Figure  40.  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 

E  =  3.0  MV/m;  a.  Variations  as  Indicated 
c  3 


2.00E+3 


J.60E+3 


1.20E+3 


8.00E+2 


4.00E+2 


Figure  41 


30  60  90  120 

time  (ns) 

Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 


E  =3.5  MV/m;  a.  Variations  as  Indicated 


d.0 


iz 


clear  time  limits— 


30 


60  90 

time  (ns) 


120 


Figure  44  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 

a .  =  0.1  V  1s  1;  E  Variations  as  Indicated 
J  C 


Figure  45.  Comparison  of  Measured  ( - )  and  Townsend 

Model  ( - )  Results: 
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Figure  4 


b.  Discussion  of  Results 


The  following  is  a  discussion  of  the  numerical  studies 
made  for  the  Townsend  corona  model: 


Variables: 


E  E-field  value  at  which  breakdown  occurs 
c 

a j  "Townsend  ionization  coefficient",  essentially  a 

measure  of  the  efficiency  of  the  avalanche  process  in 
generating  successive  electron-ion  pairs  (i.e.  pairs 
that  do  not  readily  recombine,  since  those  would  not 
be  part  of  the  corona  formation) 


^.imitations: 

Peek's  law  constrains  2.5  x  106  V/m  <  Ec  <  4.5  x  106  V/m. 

In  Townsend's  original  model,  would  equal  6.03  x 

-4  -1  -1 

10  Vs  for  this  test  geometry.  The  transient  nature  of 
the  driving  E-field  used  here  invalidates  some  of  the 
assumptions  leading  to  this  value,  and  a  dimensionless 
parameter  a.  (see  Reference  [7])  is  introduced  to  fit  the 

j  c 

model  to  the  results  obtained  experimentally.  In  the  original 

Townsend  model,  a.  =  1.0  .  Reference  [7]  makes  no 

U  c  u 
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restrictions  on  the  value  of  ajc  *  Instead,  nominal  values 
■  .u  suggested  for  . 

Variations  tested: 

E  :  1.0  X  106,  1.5  x  106,  2.0  X  106,  2.5  X  106, 

o 

3.0  x  106,  3.5  x  106  (all  in  V/m) 
a j :  lo'3,  lo“?,  10_1,  1,  5,  10,  20,  100  (all  in  V-1s-1) 

Results: 

From  the  graphs  in  Figures  36-41  ,  one  can  see  that  for  a 

constant  value  of  E  ,  larger  values  of  a  .  lead  to  greater 

suppression  of  the  induced  current  seen  on  the  line.  In  fact, 
-l  -l 

aA  -  100  Vs  provides  the  best  fit  of  the  model  to  the 

J 

general  form  of  the  experimental  data. 

Holding  a j  constant  while  varying  Ec  is  demonstrated 

in  Figures  39-45..  For  smaller  values  of  ,  variations  in 

E  show  no  effect  (see  Figures  42  and  43  ).  For  a.  >  0.1 
c  3 

_  1  _  i 

V  s  ,  som  "fanning"  of  the  response  curves  can  be  seen, 
with  the  results  for  higher  Ec  values  demonstrating  smaller 
current  suppression.  This  might  be  expected  since  as  E 

C 

increases,  the  air  will  support  greater  field  strengths  and, 
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hence,  greater  induced  currents  on  the  test  line,  before 

corona  formation  begins.  The  best  fit  of  the  model  to  the 

-1  -1 

data  may  be  seen  in  Figure  48  where  a ^  =  100.0  V  s  .  Here, 
the  range  of  E  was  extended  beyond  the  range  indicated  above 
and  a  surprisingly  good  fit  was  found  when  E  =  4.5  MV/m. 

D.  The  Mo  Corona  Model 


1.  Theory 

The  corona  model  developed  by  T.C.  Mo  and  described  in 
reference  [16]  is  different  from  the  others  in  that  the 
details  of  the  physics  of  the  corona  formation  are  used  in 
developing  the  solution.  (Only  the  salient  features  of  this 
model  will  be  summarized  here,  and  the  reader  is  referred  to 
the  original  report  for  a  more  detailed  derivation  of  the 
theory . ) 

The  Mo  model  assumes  that  transmission  line  theory  is 
sufficient  for  describing  the  propagation  of  the  current  and 
charge  along  tne  line.  Tho  telegrapher's  equations  used  by  Mo 
can  be  put  into  our  standard  form  as 


S +  £  -  Eino 


(19a) 
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and 


ai 

dz 


+  C' 


ft  -  F<v> 


(19b) 


where  L'  is  defined  as  in  equations  (2)  and  (3)  and  C'  is  a 
constant  given  by  equation  (4). 


Note  that  in  Mo’s  original  formulation,  the  factor  of  2 
in  the  logarithm  term  was  omitted.  In  addition,  in  Mo's 
development,  the  term  v  is  frequently  referred  to  as  the 
line  voltage  which,  according  to  the  discussions  in  reference 
[6],  must  have  an  additional  contribution  from  the  vertical 
component  of  the  incident  electric  field.  Actually,  the  term 
v  used  here  is  a  measure  of  the  local  charge  on  the  wire,  as 
previously  discussed. 


In  the  Mo  model,  the  effective  source  term  in  equation 
(19b)  is  shown  to  be: 


F (v)  =  -g(v)  v 


(20) 


where  g(v)  is  a  nonlinear  conductance  given  by  the  relation 


g(v) 


s 

(ac(v)  '  aw) 


In  (ac(v)/aw) 
.  In  (2h  7%) 


(21) 
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This  term  could  have  been  moved  to  the  left  hand  side  of 
equation  (19b)  and  treated  as  part  of  the  unknown,  but  it  is 
permissible  to  consider  it  as  a  ''source"  term  and  use  the  same 
calculational  algorithm  as  used  in  the  solution  of  the 
Engheta-Lee  models. 


In  equation  (21) ,  a  represents  the  effective  corona 

U 

radius  and  is  determined  in  the  same  manner  as  in  the  Baum 
model ,  namely 


a  = 
c 


Q' 


2v  a  E„ 
o  c 


2v  E 


(6) 


The  term  s  in  equation  (21)  is  an  effective 
conductivity  of  the  corona,  and  is  itself  a  nonlinear,  time 
dependent  function  of  the  local  electric  field  in  the  corona. 
In  this  model,  the  behavior  of  the  parameter  s  is  determined 
by  a  rate  equation  of  the  form 

=  «(E;,P)  I E  *  |  s  ,  (22) 

dt  pep 


where  E^  is  the  radial  electric  field  at  the  corona  boundary 
and  is  given  by 
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E* 

P 


(23) 


ln(ac/aw) 

ln(2h/aw)  '  (ac  -  ay) 

The  term  a ( E^,p)  is  an  approximate  empirical  formulation  of 

the  Townsend,  ionization  coefficient  given  by 

*  ( 

a ( E ' , p )  =  990.0  p  exp  -(2.72  x  104  p/E ' )  •  exp  (f)  (24) 

P  P 

with  f  =  (0.042  (2.72  X  104  p/E'  -  3 . 4 ) 2 ) )  ,  (24a) 

P 

where  p  is  the  ambient  air  pressure  in  mm-Hg,  and  pi&  is 
the  electron  mobility,  estimated  to  be  on  the  order  of  0.275 
m2/v-s  [ 24 ]  . 

The  solution  of  the  line  current  using  this  model, 
therefore,  involves  solving  equations  (19a)  and  (19b)  in 
conjunction  with  the  rate  equation  (22)  at  each  time  step.  In 
order  to  begin  the  solution,  it  is  necessary  to  specify  an 
initial  value  for  the  parameter  s  (denoted  by  sq) .  This  is 
a  quantity  that  may  be  specified  independently  of  other 
parameters,  and  thus,  can  be  regarded  as  a  "tuning"  parameter 
for  the  corona  model. 
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Assuming  that  the  air  pressure  is  fixed  at  621  mm-Hg 

(ambient  or  station  (i.e.  not  corrected  to  STP)  barometric 

pressure  measured  by  National  Weather  Service  at  the  tiim:  and 

for  the  altitude  at  which  the  baseline  corona  measurements 

were  made) ,  the  only  other  free  parameters  in  this  model  are 

the  line  height  h  and  the  critical  electric  field  strength 

Ec.  For  the  parametric  study  using  the  Mo  model,  the  parameter 

h  was  held  constant  at  20  meters  since  the  response  is  not 

extremely  sensitive  to  this  parameter,  and  s  and  E  were 

o  c 

varied.  Figures  49  through  60  present  the  results  predicted  by 
this  model  compared  with  those  acquired  for  the  experiment 
baseline  geometry. 
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Figure  59.  Comparison  of  Measured  ( - )  and  Mo  Corona 

Model  (—— -)  Results: 
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2.  Discussion  of  Results 


Variables: 

E„  E-field  value  at  which  breakdown  occurs 
c 

sQ  parameter  related  to  initial  value  of  air 
conductivity 

Limitations: 

Peek's  law  constrains  2.5  x  106  V/m  <  Ec  <  4.5  x  106  V/m. 

As  indicated  for  the  conductivity  model,  air  conductivity 

during  corona  is  very  difficult  to  measure  since  it  is  a 

function  of  many  parameters  including  amount  and  polarity  of 

ionization,  local  pressure,  and  relative  humidity. 

Alternatively,  it  is  a  function  of  the  mobility  and  number 

density  of  the  charge  carriers.  Here,  two  values  were  tried 

for  the  mobility  of  the  charge  carriers:  0.275  (for  electrons) 

and  0.0002  (for  ions)  [24]  (both  values  given  in 
2 

meters  /volt-second  and  representing  typical  mobilities  for 
electrons  and  ions,  respectively) .  An  empirical  formula  for 
Townsend's  ionization  coefficient  was  formulated  to  fit  the 
data  of  von  Engel  [16]. 
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Variations  tested: 


E  :  3.0  x  106,  3.5  x  106,  4.0  x  106,  4.5  x  106  (all  in 

c 

V/m) 

—  n  —A  — C 

sQ:  10  ,  10  ,  10  (all  in  mhos  or  Siemens) 

Results : 

2 

Originally,  an  ion  mobility  [24]  of  =  0.0002  m  /V-s 

was  tried  in  the  Mo  model  (see  Figures  49  -  56) .  This  value 
is  approximately  correct  for  both  positive  and  negative  ions. 
Here,  however,  it  was  anticipated  that  the  primary  charge 
carriers  were  positive.  It  was  hypothesized  that  the  test 
line  became  positively  charged  as  a  result  of  its  interaction 
with  the  simulated  EMP.  Electrons  appearing  very  near  the 
wire  were  quickly  sucked  into  the  wire  and  the  displacement 
current  noted  during  corona  is  anticipated  to  be  due  primarily 
to  positive  ions  moving  away  from  the  wire  [20].  Th^  results 
obtained  using  such  a  mobility  (see  Figures  49-56)  fit  the 
acquired  data  almort  as  well  as  did  those  of  the  conductivity 
model . 

Alternatively,  a  mobility  characteristic  of  electrons  [24] 
2 

(u  =  0.275  m  /V-s)  was  also  tried.  If  the  line  was  indeed 
© 
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charged  positively,  it  was  conjectured  that  electrons  might 
actually  be  the  charge  carriers  of  interest  [20].  In  the 
model  used,  numerical  instabilities  developed  in  the  late  time 
portions  of  the  curves  (see  Figures  57-60)  and  the  fit  to  the 
acquired  data  was  in  general  quite  poor. 

E.  The  Kudyan-Shih  Corona  Model 

1 .  Theory 

The  model  proposed  by  Kudyan  and  Shih  [13]  for 
investigating  the  propagation  of  surges  on  a  line  with  corona 
involves  modeling  the  corona  by  a  nonlinear  circuit.  The  v-i 
propagation  relations  along  the  wire  assumed  by  these  authors 
are  of  the  following  form: 

Sv  _  .  3i  (25a) 

dz  u  at 

and 

ai _ dQ'  (25b) 

ax  at 
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where  Q‘  is  the  charge  per-unit-length  on  the  wire,  and  is 
related  to  the  voltage  on  the  wire  (with  respect  to  the 
ground-plane)  as 

Q* (z,t)  =  C* (z,t)  v(z,t)  .  (25c) 

C'  is  an  effective  nonlinear  capacitance.  Kudyan  and  Shih 
postulate  that  the  circuit  model  shown  in  Figure  58  adequately 
describes  the  Q-v  relationship  along  the  wire,  and  that  a 
simple  circuit  analysis  can  be  used  to  determine  the  effective 
line  and  corona  capacitance  at  any  time  and  position. 

In  applying  this  approach  to  the  distributed  field 
excitation  of  the  line,  we  again  define  the  local  variable  v  a 
Q'/C'  and  write  the  forced  telegrapher's  equations  as  in 
equations  (9a)  and  (9b) ,  where  the  v-Q  relationship  for  the 
corona  region  is  assumed  to  be  described  by  the  circuit 
presented  in  Figure  61.  In  this  figure  there  are  two  parallel 
nonlinear  circuits  representing  the  effects  of  a  positive  and 
negative  corona.  The  batteries  labeled  as  Vc  account  for  the 
air  breakdown  electric  field  thresholds.  In  this  model,  the 
values  are  assumed  to  be  identical  for  the  positive  and 
negative  coronas.  The  resistance  and  capacitance  values  of  the 
nonlinear  elements  are  not  determined  from  first  principles, 
but  estimated  from  comparisons  with  experimental  line 
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(a)  Nonlinear  model  circuit  for  an  element  of  line  in  corona 


V(t) 

V.(t)  0  V+(t) 


(b)  At  a  given  instant,  the  circuit  will  be  in  one  of  three 

possible  states  (negative  corona,  normal,  or  positive  corona). 


Figure  61.  Circuit  Model  for  Element  of  Transmission 
Line  in  Kudyan-Shih  Corona  Model 


responses.  For  simplicity,  it  is  assumed  R^  «  R2  and 

®  C2  although  in  reality  these  values  may  be  expected  to 
differ  from  each  other. 

A3  discussed  in  [13],  the  circuit  can  be  in  one  of  three 
states,  depending  on  the  voltage  of  the  line  and  the  past 
history  of  the  line.  By  defining  a  positive  corona  model 
voltage  v+  as 

v+(t)  -  v2(t)  +  Vc  ,  (26) 

where  v2  represents  the  voltage  across  the  corona  capacitance 
C2,  and  a  corresponding  negative  corona  model  voltage 

v“(t)  »  v3(t)  -  Vc  ,  (27) 

it  is  possible  to  define  the  three  corona  states  and  the 
voltages  across  the  capacitive  elements  as  follows: 

For  v(t>  >  v+ftl :  (positive  corona  state) 

v2(t)  -  v2(t-At)  +  [v(t)-Vc-v2(t-At))  (l-e”At/(R2C2))  <28> 

v3(t)  »  v3 (t-At)  e”At/ (R3C3 ) )  .  (29) 
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For  vft)  <  v  :  (negative  corona  state) 


V2 (t)  -  v2(t-At)  e"At/(R2C2)) 


(30) 


and 

v3(t)  =  v3(t-At)  +  [v(t)-Vc-v3(t-At)]  (l-e-At^  ^R3C3^ )  .  (31) 
For  v~ ft)  <  v(t)  <  v+ft)i  (normal  state) 

v2(t)  “  V2  e"At/(R2C2))  (32) 

and 

v3  (t)  =»  v3(t-At)  e“At/ (R3C3) )  .  (33) 

With  a  knowledge  of  the  voltage  across  each  of  the 
capacitors  in  the  model,  it  is  possible  to  write  the  total 
charge  per-unit-length  stored  as 

C(t)  -  C^v(t)  +  C’v2(t)  +  C^v3(t)  .  (34) 

The  above  expression  for  the  charge  can  be  separated 
into  a  constant  capacitance  term  involving  and  the 
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remaining  terms.  This  permits  the  equatic^s  (9)  to  be 
re-written  as  though  they  contained  source  terms  as: 


fiY  .  r .  fil  .  Einc 
dz  +  h  dx  Ez 

and 


(35a) 


“  +  ci  ft  "  *  at  [  ciVt>  +  c3v3<t>  ]  <35b> 

and  the  finite-difference,  time  marching  solution  developed 
for  the  previous  models  can  be  used  to  obtain  a  solution  for 
the  line  current. 

Figures  62  through  71  present  some  of  the  results  of  a 
variation  of  parameters  study  for  this  corona  model. 
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Figure  67.  Comparison  of  Measured  ( - )  and  Kudyan/Shih 

Model  ( - )  Results: 

Ec  ■  5  MV/m;  R1#R2  as  Indicated;  C^,C2  *  10~*F 


Figure  68.  Comparison  of  Measured  i - )  and  Kudyan/Shih 

Model  ( - )  Results: 

Ec  *  5  MV/m;  R1#R2  as  Indicated;  Cj^Cj  ■  5xlO"*°F 


96 


Figure  69.  Comparison  of  Measured  ( - )  and  Kudyan/Shih 

Model  ( - )  Results: 

**  1 

Ec  «  5  MV/m;  R1#R2  as  Indicated;  C1#C2  -  10 


Figure  70.  Comparison  of  Measured  ( - )  and  Kudyan/Shih 

Model  ( - )  Results: 


«  5  MV/m;  R, ,R~  as  Indicated;  C, ,C3  - 


E 


5x1 


2.  Discussion  of  Results 


variables ; 

Ec  E- field  value  at  which  breakdown  occurs 

C1,  C2  Capacitances  in  nonlinear  model  circuit  element 

R1,  R2  Resistances  in  nonlinear  model  circuit  element 

Limitations: 

Peek's  law  constrains  2.5  x  10  V/m  <  Ec  <  4.5  x  10  V/m. 

Lacking  sufficient  information  to  make  a  realistic 
distinction  between  negative  and  positive  coronas,  it  was 
assumed  that  R^  =  R2  and  ~  C2  *  :*'s  unknown  to  wbat 

extent  this  assumption  affects  how  closely  the  model 
predictions  match  the  experimental  data. 

This  model  may  be  more  accurate  for  predicting  the  EMP 
response  of  a  conducting  line  carrying  AC,  since  the 
phenomenon  described  closely  resembles  "dielectric  hysteresis" 
as  described  in  paragraph  129,  section  2-33  of  [19],  The 
experiment  performed  utilized  a  conducting  test  line  carrying 
no  current,  AC  or  DC.  Future  tests  may  wish  to  explore  the 
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impact  of  superposi..:;  the  EMP  induced  current  on  a  current 
carrying  test  line. 

Variations...  tested: 

E  :  1.0  x  106,  2.0  x  106,  3.0  x  106,  4.0  x  106, 

c 

5.0  x  106  (all  in  V/m) 

C1  *  C2i  1.0  x  lo"9,  5.0  x  10”10,  1.0  x  10~10, 

5.0  x  lo"11,  1.0  x  10”11  (all  in  farads) 

R1  *  R2:  100'  1000»  10000,  100000  (all  in  ohms) 

Note  that  although  all  these  variations  were  tested,  the 
results  presented  are  a  subset  of  the  variations  indicated. 
Results: 

The  results  of  this  model  compared  somewhat  favorably 

with  the  experimental  data  are  presented  in  Figures  62-71 . 

For  E  <  4.0  x  106  V/m,  the  fit  of  the  model  to  the 
c 

experimental  data  is  poor  at  all  chosen  values  of  R2,  C^, 

and  C-  .  At  E_  =  4.0  ::  106  V/m,  the  fit  is  fair  for  all 

£  C 

ranges  of  and  C2  when  Ri  =  R2  *  10000  ohms  /  with  the 

best  fit  found  when  C1  =  C2  =  lO-11  farads.  At  Ec  -  5.0  x 
106  V/m,  better  fits  were  found  when  *=  R2  =  100  ohms  and 
3  x  10“10  i  C1  «  C2  i  10~9  farads  as  well  as  for  Ri  “  R2  * 
1000  ohms  and  5  x  lo"11  i  Cx  =  C2  S  10-9  farads. 
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V.  CONCLUSIONS 


Five  separate  corona  models  were  compared  against  a 
baseline  test  case  (i.e.  a  representative  example  of  the 
acquired  experimental  data)  to  determine  the  model  providing 
the  closest  fit  to  that  measured  data.  In  general,  three 
models  appeared  to  fit  the  measured  data  reasonably  well:  the 
Townsend  and  Conductivity  models  proposed  by  Engheta  et  al. 

[7]  and  the  corona  model  proposed  by  Kudyan  and  Shih  [13]. 

When  a  mobility  for  appropriate  for  ions  was  used  in  the  Mo 
model  [16],  its  predictions  fit  the  acquired  data  fairly  well 
also. 


The  Baum  model  [6]  showed  insufficient  current 
suppression  once  corona  effects  were  noted.  This  modal 
neglected  the  physics  of  corona  formation  and  assumed  that  all 
current  suppression  was  due  to  a  change  in  effective 
per-unit-length  line  capacitance  brought  on  by  the  presence  of 
corona.  The  transmission  line  was  also  assumed  to  be 
lossless.  It  is  possible  that  including  some  consideration  of 
both  the  line  losses  and  a  more  detailed  look  at  the  physics 
of  corona  and  its  formation  mighc  improve  the  predictions  of 
this  model. 
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The  model  proposed  by  Mo  demonstrated  a  remarkable  fit  to 
the  acquired  data  when  an  ion  mobility  value  was  used  for  the 
model  parameter  However,  when  a  mobility  value 

appropriate  for  electron  was  used  instead,  the  model  results 
were  unstable  and  never  approached  the  experimental  data.  This 
may  reflect  inadequacies  in  the  form  of  solution  chosen  when 
applied  to  this  particular  model  or  it  may  indicate  that  ions 
were  the  actual  charge  carriers  involved  in  corona  formation. 

The  models  demonstrating  good  fits  to  the  experimental 

data  were  the  Townsend,  Conductivity,  Kudyan/Shih  (also 

commonly  referred  to  as  the  "Hysteresis1'  model)  and  the  Mo 

model  (assuming  ions  as  primary  charge  carriers  in  the 

formation  of  corona) .  In  the  Townsend  model,  the  best  fit  is 

~1  -1 

seen  when  =  100  V  s  and  Ec  =  4.5  MV/m  (see  Figure 

48).  Recall  that  in  this  model  is,  in  essence,  a  measure 

of  the  rate  of  corona  formation.  The  Conductivity  model  most 

-5 

closely  matched  the  measured  data  when  a  =  5.0  x  10  mhos/m 

and  4.5  i  E  i  5.0  MV/m.  Here,  a  is  a  measure  of  the 

corona  conductivity.  The  Hysteresis  model  generally  fit  the 

measured  data  when  4.5  £  Ec  <  5.0  MV/m  and  =  R2 

3  4 

was  between  10  and  10  ohms.  Varying  C1  =  C2  had 

minimal  effect  within  this  range.  Perhaps  the  best  fit  of 

.  3 

this  model  to  the  acquired  data  was  for  Ri  =  R2  =  10  ohms 

—  11 

and  C^  =  C2  =  10  F.  The  Mo  model  fit  the  acquired  data 
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best  for  u  «  2  x  10“4  m2/V-s,  4.5  MV/m  <  E_  <  5.0  MV/m  and 
©  c 

s  on  the  order  of  10  mhos.  It  may  be  noted  that  none  of 
o 

these  models  fits  the  measured  data  extremely  well  in  the  time 
prior  to  corona  formation  and  immediately  after  its  effects 
are  noted. 

Clearly,  much  work  remains  to  be  done.  The  comparisons 
made  in  this  paper  indicate  a  limited  group  of  the  comparisons 
that  could  be  made.  Some  models  may  demonstrate  better  fits 
for  other  sets  of  acquired  data.  By  the  same  token,  the 
models  demonstrating  good  fit  for  the  selected  baseline  case 
may  not  continue  to  fit  other  pieces  of  acquired  data.  During 
the  experiment,  over  400  sets  of  data  were  acquired, 
representing  more  than  20  different  test  configurations  [see 
Reference  3 ] .  It  may  be  useful  to  understand  the  variations 
within  data  sets  for  the  same  configuration  and  to  check  the 
fit  of  each  model  to  an  appropriate  representative  data  set 
for  each  of  the  configurations.  These  tests  may  more  clearly 
indicate  the  strengths  and  weaknesses  of  each  model. 
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